Dramatic impact of pumping mechanism on photon entanglement in microcavity 



Alexander N. Poddubny 
loffe Physical-Technical Institute RAS, 26 Polytekhnicheskaya, 194021 St, -Petersburg, Russia 

(Dated: October 4, 2011) 

A theory of entangled photons emission from quantum dot in microcavity under continuous and 
pulsed incoherent pumping is presented. It is shown that the time-resolved two-photon correlations 
drastically depend on the pumping mechanism: the continuous pumping quenches the polarization 
entanglement and strongly suppresses photon correlation times. Analytical theory of the effect is 
presented. 
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Semiconductor quantum dots are a promising source of 
single photons and entangled photon pairs. Polarization- 
entangled photons generated during the radiative recom- 
bination of the quantum dot biexciton are now in a focus 
of intensive experimental research [IHH] ■ 

The state-of-the art approach to increase the rate of 
photon pair generation is to place the dot in the spe- 
cially designed cavity, where the frequencies of the two 
different photon modes are independently tuned to the 
biexciton and exciton resonances Both exciton and 
biexciton radiative recombinations are then increased, al- 
lowing to observe bright two-photon emission [H]. Here 
we study theoretically the quantum emission properties 
of such microcavity under incoherent pumping. We ana- 
lyze the polarization density matrix of the photon pair, 
determined by the second-order correlation function of 
the photons g^"^^ [S]. The incoherent pumping itself is an 
intrinsic feature of the light source. However, to the best 
of our knowledge, the pumping effect on the biexciton 
emission has not been theoretically analyzed yet, despite 
the extensive amount of the studies done [31 [HI [IHl [H] ■ 
Experimental design of Ref. [2] has not been addressed 
theoretically as well. 

Here we demonstrate that the entanglement is strongly 
suppressed at incoherent continuous pumping. The qual- 
itative explanation of this effect is presented on Fig. [Tj 
Fig. [ija) schematically illustrates the cascade of biexci- 
ton emission. The radiative recombination of the biex- 
citon leads to the generation of either two horizontally 
{x) polarized (red arrows) , or two vertically {y) polarized 
(blue arrows) photons in the cavity. When anisotropic ex- 
change splitting of the bright exciton state[T2] vanishes, 
these two channels have the same probability, leading to 
the completely entangled two-photon state. The situ- 
ation changes dramatically when the excitons are gen- 
erated in the quantum dot continuously. One of the 
possible mechanisms of the entanglement suppression is 
schematically illustrated on Fig. [ijb) . At the first step 
the biexciton emits x-polarized photon (dotted arrow). 
After that one ai-polarized exciton is remained in the dot. 
Due to the pumping, another exciton with y polarization 
can be generated (curved green arrow). This y-polarized 
exciton can emit photon before the z-polarized one, so 
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FIG. 1: (Color online) (a) Scheme of the two-photon emis- 
sion from a quantum dot. Red and blue arrows correspond to 
X- and y-polarized photon modes, tuned to the exciton (solid 
arrows) and biexciton (dotted arrows) resonances. The let- 
ters 0, X, Y and XX denote ground state , two bright exciton 
states and biexciton state, respectively. Panel (b) schemati- 
cally presents the breakdown of the entangled photon emission 
at continuous pumping. 

that a pair of cross-polarized photons is present in the 
cavity. Such naive analysis hints that the pumping can 
suppress the generation of entangled photon pairs and 
raises the demand for more thorough calculation. 

To elaborate on this effect we have performed a rig- 
orous simulation based on the master equation for the 
density matrix of the system [13]. We consider a zero- 
dimensional microcavity with single quantum dot. The 
Hamiltonian % can be written in the following form, 

T~L = "Hcxc + ^phot + T~L cxc— phot 7 (1) 

where the three terms are the Hamiltonians of the dot, of 
the cavity photons and of their interaction, respectively. 
The quantum dot is grown along z \\ [001] axis from 
a zinc-blende semiconductor. Only heavy-hole excitons 
are taken into account. Under the above assumptions the 
Hamiltonian of the dot reads [T^] 

■Hcxc = ^ hiL!xnx,a + fi^xxnxx ■ (2) 

a—x.y,d.d' 

The summation in Eq. ([2| is performed over x- and y- 
polarized bright heavy-hole exciton states |v'ra;_y)[12] and 
over two dark exciton states \Xd^d'), with huJcxc.a being 
the energies of these states. The singlet biexciton state 
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with the energy tiuxx is denoted as The exci- 

ton creation operators aj^ have the only non-zero ma- 
trix elements (X„|a1,|0> ^ {XX\al\Xa) = 1. Here |0) 
is the ground state of the system with no excitations 
in conduction and valence bands and no photons in the 
cavity. The operators of exciton and biexciton numbers 
Tlx, a and nxx have the only non-zero matrix elements 
{Xa\nx,a\Xa) — 1 and {XX\nxx\XX) = 1, respectively. 
The anisotropic exchange splitting of the bright exciton 
doublet is ignored for simplicity. The photon Hamilto- 
nian reads 



n 



phot 



a—x,y 



ct—x,y 



(3) 

where c\ ^ and c\ ^ are the creation operators for the 
two cavity modes with given linear polarization a = x,y. 
We neglect the polarization splitting of the modes and 
assume, that they are independently tuned to the exciton 
and biexciton resonances, hereafter we term modes 1 and 
2 as exciton and biexciton photon modes, respectively. 
This corresponds to experimental situation of Ref. [5] 
Finally, the light-exciton interaction Hamiltonian reads 



cxc— phot 



^9 (^Ic^^c, 



(4) 



a— 2:,y , 
1^=1,2 



where the interaction constant g is chosen for simplicity 
real and the same for both modes. 

The incoherent pumping leads to the generation of the 
excitons in the dot [T3]. The excitons can decay both 
radiatively and nonradiatively. The photons, created 
by exciton recombination, leave the cavity by tunneling 
through its mirrors, and are detected in experiment. To 
account for all these processes, one has to solve master 
equation for the density matrix of the system p: I13j 



dp 
dt 



(5) 



2 (2alpa„-a„aJ,p-pa„aJ,) 

(y.—x,y,d,d' 

pel 



a—x^y 
i/=l,2 



Fx 



a—x^y^d.d' 

Here the quantities Px, Tx are the exciton pumping and 
nonradiative decay rates, and Tq is the photon decay 
rate. For simplicity the polarization and spin dependence 
of these three processes is disregarded. We note, that the 
Lindblad terms in ^ describe the generation and decay 
for both exciton and biexciton states. 

In this paper we restrict ourselves to the weak cou- 
pling regime for both exciton and biexciton resonances, 
i.e. g <C Fc. We also note, that in typical experiments 



Fx <Si Fc [T51 - fT7] . Our goal is to determine the two- 
photon density matrix 

Xal3{t,T) = Ci^a{t + T)c2^l3{t) , (6) 

describing the correlations between the biexciton photon 
emitted at the time t and the exciton photon emitted at 
the time t-\- t. The angular brackets denote both statis- 
tical and quantum mechanical averaging, the constant M 
in Eq. (|6| is determined from the normalization condition 
Trp(2)^l. 

Depending on the experimental conditions, two qual- 
itatively different situations can be realized: pulsed 
pumping and continuous pumping. The procedure to de- 
termine p^^-* is different in these two cases. 

(i) Pulsed pumping. In this regime we assume that 
short single pumping pulse creates the population of exci- 
tons and biexcitons in the quantum dot. After the pump 
switches off, the excitons start to recombine radiatively. 
Assuming that the pulse duration (5 ps in experiment 
of Ref. 6} is longer then the typical energy relaxation 
times of the carriers (being on subpicosecond scale|18)). 
but still shorter than then excitonic radiative lifetime 
l/Frad = rc/(45^), being on the order of 100 ps [5], we 
can separate the calculation into two steps. First, we find 
the density matrix poj generated by the pump pulse, from 
the equation £oPo = , where the Liouvillian £o differs 
from the Liouvillian C in Eq. ([5]) by neglecting the cou- 
pling term, g = 0. Second, we consider the spontaneous 
decay after the pump is switched off. This process is 
described by another Liouvillian, £dccay, where exciton- 
photon coupling is retained but the pumping Px is set to 
zero. Thus, we obtain a spontaneous emission problem, 
with initial conditions determined by the pump . The 
two-time correlator p^'^\t^T) then formally reads [T^ 



A/'Tr 



^2/3'''2/3'= 



exp 



(7) 



(ii) Continuous pumping. In this case it assumed, that 
the excitons are continuously generated in the quantum 
dot. The balance between the exciton generation and 
decay leads to the formation of the stationary density 
matrix p, found from the stationary solution of Eq. ([5|. 
This density matrix allows us to determine stationary 
particle numbers. The two-photon density matrix ^ 
depends only on the delay r and is given, similarly to 
Eq. 0, by 



(2) 



(r) = AATr 



(8) 



The density matrix equations were solved numerically 
expanding the density matrix over the complete basis of 
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FIG. 2: (Color online) Time dependence of concurrence 
(dashed lines) and fidelity (solid lines) at continuous and 
pulsed pumping. Calculation was performed at the follow- 
ing set of parameters: kg = 30 fJ,eV, KTx ~ 15 /leV, 
hFc = 300 ^eV, hwxx - 2hwx = -5 meV, Px = O.irc. 



the quantum dot and the cavity states. The resulting ma- 
trices of Liouvillians are sparse, which makes the calcula- 
tion of the matrix exponentials feasible. Some analytical 
results in stationary regime will be presented below. 

The obtained two-photon polarization density matrix 
has the following general structure: 



1 



2(iV|| +iVj_) 



/A^ll No\ 

7V_L 

7V_L 

\No iVii/ 



(9) 



where the order of matrix elements is xiX2, Xiy2, yiX2, 
2/12/2- Here the quantities A'^y and N± determine the prob- 
ability of the detection of excitonic and biexcitonic pho- 
tons with matching and different linear polarizations, re- 
spectively. The matrix element describes the correla- 
tions of clockwise and counter-clockwise circularly polar- 
ized photons. The structure of the density matrix Q is 
similar to that obtained in Ref. However, in our case 
the suppression of the entanglement is solely related to 
the pumping and the quantity N(j is real. The concur- 
rence of the state (|9|, quantifying the entanglement [3], 
reads 



C 



max 



Nil - 



(10) 



The results of calculation are presented on Fig. [2] For 
pulsed pumping one has Nfj = N^^ = 1 and iVj^ — 
0, which corresponds to the completely entangled Bell 
state, [5] where fidelity and concurrence are both equal to 
unity. For continuous pumping both concurrence and fi- 
delity are smaller then unity, and are suppressed at large 
delays r. After a certain delay r* the concurrence is 
zero, i.e. the state is not entangled. Entanglement sup- 
pression is directly related to the incoherent nature of 
the pumping. Indeed, excitons are constantly generated 
in the dot, and then eventually emit photons. As soon 
as the delay becomes larger, than the dot repopulation 
time, the temporal coherence of the photons is lost. At 
very large delay one has iVy = N± = 1/4 and Nfj = 0, 
which corresponds to completely independent photons. 
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FIG. 3: (Color online) Stationary particle numbers and cor- 
relators as function of pumping, (a) The total numbers of 
bright excitons (solid curves) and biexciton (dotted curves) 
photons as functions of pumping, (b) The numbers of exci- 
ton (solid curves) and biexciton (dashed curves) photons as 
functions of pumping, (c) Concurrence (solid curves) and fi- 
delity (dotted curves) of the two-photon pair at zero delay, 
r = 0. (d) Concurrence lifetime t* , divided by exciton ra- 
diative lifetime l/Frad- Thick curves are results of numerical 
calculation, thin curves present the analytical results of Sup- 
plemental Materials. Calculated at the same parameters as 
Fig.i 



Fig.[3]presents more detailed analysis of stationary cor- 
relators. Interestingly, even the stationary entanglement, 
calculated at r = 0, is suppressed by the pumping. First 
two panels present the calculated bright exciton, biex- 
citon and photon occupation numbers as functions of 
pumping. Thick lines correspond to the results of the 
numerical calculation. Thin lines present the analytical 
results in weak coupling regime, obtained in Supplemen- 
tal Materials. Here we briefly summarize them. Neglect- 
ing the exciton-photon coupling one can get stationary 
exciton population numbers 



Px 



APx + Pi/r 



Nxx - ^Nx,x, 
1 X 



(11) 



where a = x,y,d,d' . Eq. (Ill indicates that at small 



pumpings Nx increases linearly with pumping, and iVxx 
increases quadratically. At high pumping the dot is in 
biexciton state, i.e. A^xx 1 and A^x 0. The photon 
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(12). 

numbers A^^ ^ in the weak coupling regime are given by 



where 



(1) 



27rg^Nx,o.gu Tc^g^ = 2^g2^xxe2 (12) 



9i 



A^c + Tx + 5Px) ' 



92 



it{Tc + srx + Px) ' 



and a — x,y . The right hand sides in Eqs. (12) are 



the rates of photon generation due to exciton recombi- 
nation. They are presented in the form similar to the 
Fermi Golden rule and have simple physical meaning: 
the rate is proportional to the population of the corre- 
spondent state of the dot, to the square of the inter- 
action constant g and to the effective density of states 
g. The latter is quenched at high pumping due to the 
fermionic nature of the electrons, which typically leads 
to the pumping- induced linewidth broadening [121 US]- 
Thus, at small pumping the exciton and biexciton pho- 
ton numbers are proportional to the population of the 
corresponding quantum dot states, and at high pumpings 
they are suppressed by the pumping-induced dephasing. 
This explains the behavior of the curves on Fig. [Sja. 

At the vanishing pumping one has N(j{t = 0) = 
A^ll(r = 0) = 1 and Nj_{t ^ 0) = in Eq. which 
corresponds to the completely entangled Bell state with 
concurrence C = 1, see Eq. ([To]). The pumping leads to 
the growth of N± and to the suppression of A^q. This 
quenches the entanglement, in agreement with Fig. [sj^c). 
Fig. |3jd) demonstrates that the concurrence lifetime r* , 
i.e. the time, during which the state remains entangled, 
decreases with pumping, because the dot is faster re- 
populated. We note, that the high sensitivity of two- 
photon correlations to the pumping is a rather general 
effect, known, for instance, in superradiant emission of 
the cavities with several resonant quantum dots [21] or 
for single-dot lasers [22]. The main result of Fig. |3] IS 
that the entanglement degree is substantially quenched 
when the pumping is large enough to make the numbers 
of exciton and biexciton photons comparable. 

To summarize, we have put forward a general the- 
ory of the entangled photon generation from the zero- 
dimensional microcavity with embedded quantum dot 
under incoherent pumping in the weak coupling regime. 
We have demonstrated that the time-dependent polariza- 
tion density matrix of the entangled photon pair is very 
sensitive to the mechanism of pumping, i.e. pulsed one 
or continuous one. Analytical theory of this effect has 
been presented. Important goal of the future studies is 
to analyze the strong coupling regime, up to now realized 
only for the exciton resonance [13 [TS] . 
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Supplemental material 
SI. Stationary exciton and photon numbers 

Now we obtain the stationary occupation numbers of 
the quantum dot states and the stationary photon num- 



Assuming, that the binding energy of the biexciton 
ioxx — 2a;x is much larger than Px, Tc, Fx and <?, we 



write down Eq. (S5 1 to find biexciton-related term in pi: 



1 



(Fc + 5rx + Px)pi 



ig^xx 



6^|0)(0|ac2,a; + /i.c. (S9) 



bers, Eq. (|12|) and Eq. (|llj). We consider the weak cou- Finding pi and calculating the commutator (fSSl) again, 



pling regime, and assume that the additional condition obtain the generation rate of the biexciton photons. 



Tc > Trad, i.e. rcFx > ^ 



(SI) 



holds. Our goal is to expand the density matrix in powers 
of the coupling constant 



P = Pq + Pi + 



Pj « 



(S2) 



The lowest order contribution to exciton number is deter- 
mined by po and does not depend on 5, while the photon 
numbers are proportional to g^. To perform the expan- 
sion we express the total Liouvillian C of the system ([S]) 
as Co+Ci, where Co is the Liouvillian neglecting exciton- 
photon interaction, and 



Cip — ^ [p, "Hcxc-phot] 



(S3) 



The zero-order contribution to the stationary density ma- 
trix is found from the equation 



CqPq = . 



(S4) 



Each following term of the expansion ( S2 ) is determined 
by the recurrence relation 



^OPj — ^[Pj-l, 'Hcxc-phot] 



(S5) 



Let us first solve Eq. (S4). Obviously, the matrix po is 



diagonal, with the only non-zero matrix elements 

{X,a\po\X,a)^Nx,o. = Nx, (S6) 
{XX\po\XX) = Nxx , 

(0|po|0) - 1 - TVxx - 4iVx , 

where we took into account the normalization condition 
Trpo = 1. Substituting the density matrix in the form 



( S6 ) into Eq. ( S4 ) we obtain the following system of linear 



equations: 

- (r^ + Px)Nx + TxNxx + Px{l - Nxx - 4A^x) - , 

- r^^xx + PxxNx = . (S7) 



Solution of this system yields Eq. ( 11 ) 



Now we proceed to the calculation of the stationary 
photon numbers. We consider the case of the x-polarized 
biexciton photons as an example. Calculating the com- 



mutator Eq. (S5) for the term Nxx\XX){XX\, entering 
Po, we get 



-[po,K 



exc— phot J 



igN: 



XX 



6t|0)(0| 



ac2 . 



h.c. 



(S8) 



standing in the r.h.s. of Eq. (12). We see, that Eq. (S9) 



yields the density of states g2- The value of p2 is in- 



versely proportional to the decay rate in l.h.s. of Eq. ( S9 1 



Eq. (12 1 is in fact the particular case of Eq. (S5) to find 



the matrix p2, averaged over the states of the dot. Anal- 
ogous procedure yields the generation rate of the exciton 
photons. The consideration differs only in the decay rate 
of the intermediate state of type a^|0)(0|ci, equal to 7r/gi. 



S2. Stationary two-photon density matrix 

In this Section we present the explicit results for the 
two-photon stationary density matrix The calcula- 
tion procedure is generally the same as that presented 
above to obtain the stationary photon numbers. How- 
ever, it requires expansion of the stationary density ma- 



trix up to the fourth term in the series over g (S2|, pro- 
portional to g'^. The calculation is therefore much more 
tedious, but still feasible. Thus, we will only present the 
results: 



iVii 



(SIO) 



[ANi,xx + A'N2.\\+BNxx] , 



N± 

where 

A = 



'^2,x'A,y'^l,y'^2,x 



T^g 



[A'N2,c 



[AX 



1,XX 



BNxx] , 
(Sll) 

- A'N2,±] 

(S12) 



^(SFc + SEx + Px) 
ng'^AA' 



A' 



7r(3rc 



G 



2r> 



2P^ 



Tr.N, 



c^'c.x 



5Px) ' 
(S13) 

(2) 



20 



Px 



The correlator Ni xx is readily found from the linear 
system 

-(Ec + 4Px)iViG + 4rxiViG - -27:g^Nx,xgi, (S14) 
-(Fc + Fx + Px)A^ix + iViG^x + NixxTx = 0, 
-(Fc + 4rx)Nixx + ^NixPx = . 
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Here the quantities Ni^q, iVi^x and Ni^xx determine the 
probabihties to find one exciton photon and the dot in 
ground, excitonic and biexciton states, respectively. The 
correlators A^2,|| and N2^± are found from 

-{Tc + 4Px)iV2,G + rx(iV2.|| + N2.± + 2Ar2dark) = 0, 

(S15) 

-(Fc + Fx + Px)N2\\ + N2^gPx + A^2,xxrx = 

- 2ng'^Nxxg2, 
-{Tc + Tx + Px)N2± + N2.gPx + N2^xxrx = 0, 
-(Fc + Fx + Px)N2M + N2^gPx + N2,xxrx = 
-{Tc + 4rx)iV2,xx + (iV2,|| + N2^i_ + 2N2,d)Px = . 

Here N2,g, ^2,|l, ^2,j_, ^^2,d, and iV2,xx are the proba- 
bilities to find one x-polarized biexciton photon and the 
dot in the ground state, x exciton state, y exciton state, 
dark exciton state and biexciton state, respectively. 
Let us comment on the calculation details. The equa- 



the denominators given by the pumping-dependent de- 



cay rates, see Eq. (S13). Final equation to find the two- 
photon density matrix corresponding to j = 4, is 
trivial, since we are interested only in the value of p^^-* 
traced over the states of the dot. This equation yields the 
prefactor cx 1/Fc in Eqs. (S10)-(S12|. However, equa- 



tions (S14) and (S15) to find the components of the den- 
sity matrix p2 become nontrivial. These equations ex- 
plicitly depend on the pumping. For instance, the quan- 
tity N2^± describes the probability, that the j/-polarized 
exciton is generated in the quantum dot, after the x- 
polarized biexciton is emitted. Such process contributes 
to the correlator ( S12| and destroys the entanglement. 
It is illustrated on Fig. [TJd. Another similar pumping- 
induced process is the the absorption of the biexciton in 
the dot after the emission of the exciton photon. It is 
described by the correlator A^i,xx and also contributes 
to Eq. (S12) . At the vanishing pumping both these pro- 



tions for pi and are similar to Eq. (S9) and yield 



cesses are quenched, which leads to the totally entangled 
state with N« = N(j, Nj^ = 0. 



